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ABSTRACT 

Based on the Millennium Simulation we examine assembly bias for the halo properties : shape, 
triaxiality, concentration, spin, shape of the velocity ellipsoid and velocity anisotropy. For consistency 
we determine all these properties using the same set of particles, namely all gravitationally self-bound 
particles belonging to the most massive sub-structure of a given friends-of-friends halo. We confirm 
that near-spherical and high-spin halos show enhanced clustering. The opposite is true for strongly 
aspherical and low-spin halos. Further, below the typical collapse mass, M*, more concentrated halos 
show stronger clustering whereas less concentrated halos are less clustered which is reversed for masses 
above M*. Going beyond earlier work we show that: (1) oblate halos are more strongly clustered than 
prolate ones; (2) the dependence of clustering on the shape of the velocity ellipsoid coincides with 
that of the real-space shape, although the signal is stronger; (3) halos with weak velocity anisotropy 
are more clustered, whereas radially anisotropic halos are more weakly clustered; (4) for all highly 
clustered subsets we find systematically less radially biased velocity anisotropy profiles. These findings 
indicate that the velocity structure of halos is tightly correlated with environment. 

Subject headings: cosmology: theory - cosmology : dark matter - methods: N-body simulations 



1. INTRODUCTION 



The dependence of halo clustering on a second pa- 
rameter in addition to mass is now generally referred 
to as assembly bias. Analytic models predict t he clus- 
tering of halos to depend on their mass alo ne (jKaiserl 
198llCole fc KaiseH[l98llMo fc Whit3[T99l . However, 



Gao et al. I (I2005D and various subsequent studies showed 



that clustering also depends on other halo properties 
like formation ti me, concentration, substructure content , 
spin and shape (Harker et al.ll2006HWechsler et al.l 



Bett et all l2007t iGao fc White! l2007l: iJing et all IZUUYI: 
Maccio et al.l 120071 : IWetzel et al.l 120071 : lAngulo et all 
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It soon became clear that the dependency of clustering 
on the various parameters does not simply follow from 
the relation among those parameters. Here are two ex- 
amples: (1) Earlier formed halos have higher concen- 
trations (e g..lNavarro et al.|[T997t IWechsler et"aI1 l2002l : 
iNeto et al.l l2007t ) nevertheless the dependence of clus- 
tering on formati on time and conc entration show differ- 
ent behavi o r (cf., iJing et al]|2007f ). (2) As discussed in 
iBett et all ()2007D there is a weak correlation between 
spin and shape, more spherical halos show on average 
a slightly lower spin parameter. However, the bias be- 
havior is opposite to what one would naively derive from 
this correlation. The most nearly spherical halos and 
the halos with highest spins a re clustered above average 
(see also, iGao fc White! l2007l ). These examples suggest 
that the ranking according to any given parameter, like 
concentration, spin or shape cannot be converted in any 
simple way to the approximate ran king according to an - 
other parameter. In this context !Croton et al.! (j2007! ) 
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speculated (but did not demonstrate) that there may be 
a more fundamental parameter capable of uniquely pre- 
dicting halo clustering. 

Simple implementations of extended Press-Schechter 
and excursion set models (!Press fc Schechter! !l974! : 
!Bond et al.! !l99l!: !Lacev fc Cold !l993l ) commonly used 
to estimate halo statistics do not predict a dependence 
of formation time on environment. However, using a 
mass filter in c onfiguration space rather than in k-space 
!Zentneil (|2007l ) demonstrated that halos in denser envi- 
ronments do form later independent of halo mass. At 
the high mass en d this agrees with find i ngs from N-body 
simulations fe.g.. !Wechsler et al.![200^ : !Jing et al.!!2007f) 
but is op posite to the beh avior observed at the low mass 
e nd (e.g.,!G ao et al!!200j. 

!Dalal et a l. (20081) argued that for high halo masses as- 
sembly bias is expected from the statistics of the peaks of 
Gaussian random fluctuations and at low masses it arises 
from a subpopulation of low-mass halos whose mass ac- 
cretion has ceased. It is unclear whether whether other 
aspects of assembly bias can be explained in this way. 

Several other studies have investigated the depen- 
dence of halo formation 

rates on environment [e 

200l iSheth fc TormenI l200 



times or , similarly, merger 



!Gottl5ber et all !2001 



iFakhouri fc Mai l2009bl la 



Hahn et al.! !2009f )" Although, slightly different density 



estimators, like over density in a sphere or mark correla- 
tion functions, are employed it is generally agreed upon 
that halos less massive than ~ 10^'^ /i~^M0 which reside 
in high density regions form earlier compared to those 
with equal mass but located in less dens regions. 

Most of the previous studies examined the impact of 
the assembly bias on galaxy clustering statistics. This 
is not the prime focus of this work. We are mainly con- 
cerned about the interplay between clustering and the 
dynamical structure of halos. To describe the dynam- 
ical structure of halos we use a variety of parameters, 
namely: the mean shape, the mean triaxiality, the spin, 
the non-dimensionlized total velocity dispersion (a mca- 
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sure of concentration), the shape of the global velocity 
dispersion tensor, and the mean radial/tangential veloc- 
ity anisotropy. We also explore whether subsets of ha- 
los which show assembly bias share any common prop- 
erty. The most obvious we find is the behavior of velocity 
anisotropy. On average, all highly clustered subsets ap- 
pear to have less radially biased internal motions. The 
opposite holds for less clustered subsets. 

2. METHODOLOGY 

2.1. Simulation 



The Millennium Simulation (|Springel et all 120051) 

adopted concordance values for the parameters of a 
flat A cold dark matter (ACDM) cosmological model, 
f^dm = 0.205 and ilb = 0.045 for the current densities in 
CDM and baryons, h — 0.73 for the present dimension- 
less value of the Hubble constant, erg = 0.9 for the rms 
linear mass fluctuation in a sphere of radius 8 /i~^Mpc 
extrapolated to z = 0, and n = 1 for the slope of the pri- 
mordial fluctuation spectrum. The simulation followed 
2160'^ dark matter particles from z = 127 to the present 
day within a cubic region 500 ft,~"'^Mpc on a side resulting 
in individual particle masses of 8.6 x 10^ h~^MQ. The 
gravitational force had a Plummer- equivalent comoving 
softeni ng of 5 /t~|kpc. The Trce-PM N-body code GAD- 
GET2 (iSprineell [200l was used to carry out the sim- 
ulation and the full data were stored 64 times spaced 
approximately equally in the logarithm of the expansion 
factor. 

2.2. Halo sample and properties 

The halos are found by a two-step procedure. In the 
first step all collapsed halos with at least 20 particles are 
identified using a friends-of-friends (FoF) group-finder 
with linking parameter b = 0.2. These objects will be 
referred to as FoF-halos. Then pos t-processing with the 
substructure algorithm SUBFIND (jSpringel et al.ll200lh 
subdivides each FoF-halo into a set of self-bound sub- 
halos. Here we only consider the most massive sub-halo 
within each FoF-halo which on average comprises 80% 
of the total mass of the FoF-halo. We refer to these sub- 
halos simply as halos. As is common practice we chose 
the position of the most bound particle of that sub-halo 
as center of the halo. All halo properties are consistently 
computed based on particles assigned to that (most mas- 
sive sub-) halo. That way possible confusion due prop- 
erties which are computed based on different halo def- 
initions can be excluded. For instance, commonly the 
spin parameter is determined withn a spherical volume 
whereas the shape is based on all particles associated 
with the FoF-halo. In this study, however, both quanti- 
ties are computed based on the same set of particles. 

As is standard, the shape, s = a/c^ the triaxiality, 
T = (a^ — 6^)/(a^ — c^), and the shape of the veloc- 
ity ellipsoids, Sv — iv/cv arc computed based on the 
eigenvalues (a(v) > 6(v) > "-(v)) of the second moments 
of the spatial and velocity distributions of all halo par- 
ticles. The velocity anisotropy parameter is given by 
/? = 1 — 0.5CTj/cr^, where and at are the radial and 
tangential velocity dispersion based on all halo particles. 
Since the conventional determination of the spin parame- 
ter and the concentration is based on the particle content 
within a sphere and not on self-bound particles we intro- 
duce slightly different quantities to replace them. The 



spin, I, is computed as the absolute value of the total an- 
gular momentum multiplied by the Hubble constant H 
to the one third power and divided by the total mass M 
to the p ower 5/3. This is clo sely related to the formula 
given m iBullock et cdl (|200lh . As quantity measuring 
concentration, we use C = a/{HMy/^ where a is the to- 
tal velocity dispersion of the halo. Concentration and a 
are positively correlated, i.e. high a values correspond to 
high concentrations and vice versa. The powers of mass 
in the denominators of I and C are chosen to compensate 
for the intrinsic mass dependence of these quantities. 

2.3. Determining the bias factor 

The determinat ion of the bia s facto r, b, follows the 
approach given in iGao fc Whiti (|2007t l. It is computed 
as the relative normalization factor which minimizes the 
mean square of the difference log — log 6^„im for four 
equal width bins in log r spanning the comoving separa- 
tion range 6 < r < 20 h~^Mpc. Here, ^hm is the halo 
mass cross-correlation function and ^mm is the mass au- 
tocorrelation function. Compared to estimators based 
on the halo autocorrelation function this estimator has 
highly improved noise characteristics because of the large 
number of dark matter particle s available. As sugg ested 
by standard halo bias i nodels (|Mo fc White! 19961 ) and 
verified numerically by iGao fc White! ( 20011 ) this ap- 
proach is equivalent to estimators based on halo auto- 
correlation functions. The large-scale bias depends on 
mass and redshift through the equivalent peak height 
i^(M, z) = Sc{z)/a{M, z), where cr(M, z) is the rms linear 
overdensity within a sphere which contains the mass M 
in the mean, and Sc{z) is the linear overdensity threshold 
for collapse at redshift z. The typical collapse mass, M*, 
corresponds to = 1. We will use the equivalent peak 
height to seamlessly combine results derived at redshifts 
z = 0, 1,2 and 3. 

3. RESULTS 

3.1. Assembly bias based on consistent halo properties 

Fig.[T]shows the assembly bias based on halo properties 
as defined in $ 12.21 Symbols (squares, diamonds, trian- 
gles) display bias factors for all halos. Solid and dashed 
lines give the bias based on the upper and lower 20 % tiles 
of the distributions of each specific property as indicated 
by labels in the lower right corners. We use solid red 
lines for the more strongly biased tail and dashed blue 
lines for the more weakly biased tail. These line styles 
are adopted in Fig [2] to ease comparison. Given that 
C and I correspond closely to conventional concentra- 
tion and spin we confirm the results reported in previous 
work. However, with triaxiality, T, shape of the velocity 
ellipsoid, Sv, and the velocity anisotropy, /3, we add three 
more properties which show significant assembly bias. 

The displayed dep endency of bias o n shape, s, repro- 
duces the results bv iBett et al.l ()2007l ). For all masses, 
more nearly spherical halos are clustered above, and 
strongly aspherical halos below average. The behavior 
of the triaxiality, T, indicates that prolate halos are clus- 
tered below average and oblate halos are more strongly 
clustered than average. The results for the scaled veloc- 
ity dispersion, C, which here replaces the conventional 
concentration, indicate that for i' ^ I high concentra- 
tion halos arc more and low concentration halos are less 
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Fig. 1. — Bias factor as a function of halo mass and halo properties. Halo mass is given parametrically through the equivalent peak height 
^ = 5c/cr. The additional properties in the six panels are; shape (s), Triaxiality (T) , concentration (C), shape of the velocity ellipsoid 
(sv), the spin (L) and the velocity anisotropy (/3). Symbols are bias factors defined for all haloes in a given mass bin. The various shapes 
(squares, diamonds, two kinds of triangles) indicate the redshifts (0,1,2,3) of the snapshots used to determine them. The same redshifts 
are used for the computation of the solid curves which arc bias factors for haloes in the lower and upper 20 % tails of the distributions of 
the particular property concerned. 



clustered. This behavior reverses sign f or > 1. Th e 
same trends have been reported u r Woclis ler et al.l (I2006D. 
iGao fc Whitl ((20071 ). [Jhig et al.l(l2007 l) and lWetzet'eTaLl 
(I2007fl. Th e par ticularly good agreement with Fig. 2 in 
iJing et al] (|2007| ) based on the conventional concentra- 
tion increases confidence that C can represent concen- 
tration in the context of assembly bias studies. The 
shape of the velocity ellipsoid, Sv, behaves similarly to 
the shape of the spatial distribution, s. The dependence 
of clustering on spin, I, is comparable to the results in 
iBett et all ll200l and IGao fc Whiti (l2007l ). Finally, ha- 
los with low velocity anisotropy, (3, are more clustered 
and the opposite holds for halos with strongly radially 
biased velocities. The signal based on f) is the most pro- 
nounced among all shown here indicating that velocity 
related quantities arc tightly correlated with the cluster- 
ing strength. 

We also obtained results for the dependence of clus- 
tering on f ormat ion time (not shown) as examined in 
IGao et aTl (|2005[ ). As formation time for their f of -halos 
they used the time when half the final was acquired. For 
a variety of reasons we here restrict ourselves to self- 
bound main subhalos which are simpler and better de- 
fined objects. Th us we computed t he formation times 
in the same way as iGao et all ()2005() except that we use 
main subhalos. It turned out that the effect is qualita- 
tively similar to that found earlier but is much weaker. 
We find only a ~ 10 per cent increase of the bias fac- 



tor for early formed halos at masses corresponding to 
1/ < 1. Different factors like the orbital parameters or 
the concentrations of the merging systems imply a cer- 
tain stochasticity of the sub-halo masses which propa- 
gates into the determination of th e format i on tim es and 
reduces the assembly bias signal. iLi et al.l (|2008f ) noted 
that the definition of the formation time has a substantial 
impact on the strength of the assembly bias, confirming 
that the change from fof- to sub-halo accretion histories 
is responsible for the reduction of the assembly bias sig- 
nal. 

3.2. Do highly clustered halo subsets share similar 
properties ? 

To find out whether highly clustered halo subsets share 
similar properties we compute stacked profiles of the ha- 
los belonging the 20 % tail subsets of a given property, 
and compare them with the stacked profiles of all halos. 
We examine the difference in density, velocity disper- 
sion, phase space density and velocity anisotropy profiles. 
Density and velocity dispersion profiles show no unique 
trend with the bias behavior. 

The average phase space density profiles determined 
by p/a^ indicate that more strongly clustered subsets, 
solid red lines in Fig. [Tl tend to lie above the average 
profiles based on all halos in a given mass bin. However, 
for some properties, like spin and velocity anisotropy, the 
differences are small and dependent on radius. 
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Fig. 2. — Difference between the average velocity anisotropy profiles, A/3, of halos in the 20 per cent tails of the according distribution 
and all halos at a given mass. The same halo properties are used as in Fig. [T] From top to bottom they are: shape, triaxiality, scaled 
velocity dispersion, shape of velocity ellipsoid, spin and velocity anisotropy. Columns correspond to every second symbol for < 1 in 
Fig. [T] The two short vertical lines indicate the averages of the radii including a mean density of 200 times the critical and 200 times the 
cosmic mean density, respectively. The averages are based on all halos in a given mass bin. Averages based on halos in the tails of the 
distributions show only marginal differences. In the rightmost column the outer marks exceed the plot range. Signals beyond the marks 
are dominated by effects of the boundaries of our halos. 



IMotivated by stu dies of [Hansen et al.l (|2006[ ) and 
I Schmidt et al.l ()2008[ ) wc also probe the behavior for 
phase space densities computed as p/aj?, where de- 
notes the radial velocity dispersion. For masses larger 
than the average phase space profiles of subsets which 
show stronger clustering (solid red lines in Fig.[T]) system- 
atically lie above the average profiles of all halos. Profiles 
based on less clustered subsets fall below the average 
profiles based on all halos in a given mass bin. How- 
ever, for halo masses below the reversal of sign seen 
in concentration (vel ocity dispersion) is not reproduced. 
ISchmidt etlHI (|2008D arg ued that there is no empirical 
justification (but see also iDehnen fc McL aughlin 2005) 
for the particular definition of phase space density used 
here, and we find that changing the exponent of the ve- 
locity dispersion definitely has an impact on the correla- 
tion between clustering and phase space profiles. Due to 
such uncertainties we leave the discussion at this point. 
Nevertheless, we note there is a systematic trend for more 
clustered halos to have higher phase space densities. 

In contrast to the phase space profiles the velocity 
anisotropy profiles, (3{r) = 1 — 0.5(T^(r)/(T^ (r), are un- 
ambiguously defined (crir) and (Tt(r) denotes the radial 
and the tangential velocity dispersions averaged in spher- 
ical shells). Fig. [5] shows the difference, A/3, between the 
stacked velocity anisotropy profiles of all halos and the 
halos belonging to 20 % tails of the distribution of a given 



property. These properties are, from top to bottom: 
shape, triaxiality, scaled velocity dispersion, shape of the 
velocity ellipsoid, spin and total velocity anisotropy. The 
columns show various mass bins corresponding to every 
second mass bin with v < 1 in Fig [T] The comparison 
with Fig. [1] indicates that more clustered subsets (solid 
red lines) show on average lower (3 profiles. The opposite 
holds for the less clustered subsets (dashed blue lines). 
For the scaled velocity dispersion, C, the difference is 
strongly reduced, however the reversal of the sign as ap- 
parent in Fig. [1] is not explicitly reproduced. But apart 
from that, there is a clear correlation between clustering 
behavior and velocity anisotropy. 

One possible explanation for the more nearly isotropic 
global velocity distribution of more clustered halos is 
that the impact parameters of the merging sub-halos 
are larger due to gravitational deflections short time be- 
fore accretion. Fakhouri & Ma (2009a) have shown that 
mergers arc the dominant mode of accretion in high den- 
sity regions thus they presumably determine the velocity 
structure of these halos. On the other side in lower den- 
sity regions the gravitational field is dominated by the 
halo itself and accretion occurs in a more radial fashion 
leaving the observed imprint on the velocity structure of 
the halo. 

4. CONCLUSION 
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Taking advantage of the large volume and the supe- 
rior resolution of the Millennium Simulation we reassess 
assembly bias for various halo properties. We conclude 
with a recapitulation of the main results: 

1) Based on consistent determinations of various halo 
properties, like shape, spin and concentration, we con- 
firm the results reported in previous studies. More nearly 
spherical and high-spin halos are clustered above and 
aspherical and low-spin below average. Below the typ- 
ical collapse mass, A/*, more concentrated halos show 
stronger clustering whereas less concentrated halos are 
less clustered. This reverses for masses above M*. The 
clustering cannot simply be explained by the relations 
among the properties. 

2) To the halo properties which have already been 
shown to correlate with clustering behavior, we have 
added the triaxiality, the shape of the global velocity 
ellipsoid and the velocity anisotropy parameter. Oblate 
halos, halos with a more nearly isotropic global veloc- 
ity distribution and halos with weak velocity anisotropy 
are more clustered than average. Contrary properties 
reduce clustering. Very prominent signals are found for 
the velocity-related properties, indicating that halo ve- 



locity structure is tightly correlated with the clustering 
behavior. 

3) We also showed that independent of the halo prop- 
erty (shape, triaxiality, spin, concentration) the velocity 
anisotropy profile of the more clustered subsets lies sys- 
tematically below the average profile, whereas less clus- 
tered subsets show profiles, indicating more radially bi- 
ased internal motions. 

Our findings show that the internal velocity structure 
of halos is strongly influenced by environment. Halos are 
not "universal" and their internal properties (the devia- 
tions from universality) are related in complex ways to 
their environment. 
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